PT symmetric lattices with a local degree of freedom 
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Recently, open systems with balanced, spatially separated loss and gain have been realized and 
studied using non-Hermitian Hamiltonians that are invariant under the combined parity and time- 
reversal (VT) operations. Here, we model and investigate the effects of a local, two-state, quantum 
degree of freedom, called a pseudospin, on a one-dimensional tight-binding lattice with position- 
dependent tunneling amplitudes and a single pair of non-Hermitian, "PT-symmetric impurities. We 
show that if the resulting Hamiltonian is invariant under exchange of two pseudospin labels, the 
system can be decomposed into two uncoupled systems with tunable threshold for VT symmetry 
breaking. We discuss implications of our results to systems with specific tunneling profiles, and 
open or periodic boundary conditions. 



Introduction: A non-Hermitian Hamiltonian H ^ 
that is invariant under the combined parity and time- 
reversal (VT) operations is called VT symmetric. Since 
the groundbreaking discovery of such Hamiltonians in 
continuum models fifteen years ago [1], significant re- 
search has been carried out to identify and characterize 
the properties of PT-symmetric Hamiltonians that, typ- 
ically, are decomposed into a Hermitian kinetic term and 
a non-Hermitian, T-'T'-symmetric potential term, V{x) — 
V*{-x) ^ V^x) [HIS]- The spectrum cx of such a non- 
Hermitian Hamiltonian is purely real over a region of the 
parameter space, called the 7T-symmetric phase; in this 
region, its (non-orthonormal) eigenvectors are simul- 
taneous eigenfunctions of the VT operation. For param- 
eters outside the PT-symmetric phase, the eigenvalues of 
the Hamiltonian occur in complex conjugate pairs, and 
due to the anti-linear nature of the time-reversal operator 
T, the corresponding eigenfunctions are not simultane- 
ous eigenfunctions of the VT operation. This emergence 
of complex eigenvalues is called PT-symmetry breaking. 
PT-symmetric Hamiltonians are ideally suited to model 
non-equilibrium phenomena that transition from a quasi 
steady-state behavior (PT-symmetric phase) to loss of 
reciprocity (broken PT-symmetry) [4, 5 . 

In the past three years, experiments on coupled opti- 
cal waveguides [MS], coupled electrical circuits TU], and 
coupled pendulums [TTj have shown that instead of being 
a mathematical curiosity, PT-symmetric Hamiltonians 
represent open (quantum) systems with spatially sepa- 
rated, balanced, loss and gain. The discrete nature of 
these systems has also sparked new interest in the prop- 
erties of PT-symmetric tight-binding lattice models with 
different topologies [12 [T^] ; such lattice models are most 
readily realized in evanescently coupled optical waveg- 
uides [m dnj. Recent theoretical work has led to the 
identification of robust and fragile PT-symmetric phases 
in a lattice with open boundary conditions |16H19'. tun- 
able PT-symmetric threshold in a lattice with periodic 
boundary conditions 20 , and substantially strength- 
ened PT-symmetric phase in finite lattices with position- 
dependent tunneling profile [211 [22]. AH of this work is, 
however, restricted to systems in one spatial dimension 
where the parity operation is defined as P : a; — ^ —x in 



the continuum case (with a suitably defined origin) and 
V : k —i' k = N+l—k in a lattice with N sites. In particu- 
lar, the properties of PT-symmetric Hamiltonians in two 
(or higher) dimensions have been barely explored [53] . 

In this paper, we investigate PT-symmetric lattices 
with a local, two-state, quantum degree of freedom la- 
beled by a pseudospin a = ±1. We present a class of 
models that can be mapped onto one-dimensional lat- 
tice models that have been investigated in the past, and 
thus are solvable in a straightforward manner. Such 
degree of freedom can represent, for example, two or- 
thogonal polarizations of a mode in a single elliptical 
waveguide 24J in an array of coupled elliptical waveg- 
uides. Thus, although we use the term "pseudospin" 
to denote this degree of freedom, we emphasize that its 
time-reversal properties are unspecified. By using phys- 
ically motivated PT-symmetric potentials (at only two 
sites) and tunneling amplitude profiles, we show that the 
local degree of freedom leads to a robust, tunable PT- 
symmetric phase. Although our results are applicable to 
general PT-symmetric systems, in the following, we use 
a language that is applicable to coupled optical waveg- 
uides [TS]. 

Tight-binding Model: We consider a lattice of N coupled 
waveguides with open boundary conditions, described by 
the following tight-binding Hamiltonian, 

Ho = -Y^ [4 r(i)a,+i -I- al+iTt (^)a,] . (1) 

Here aj^ — (a^^,aj, _) are the creation operators for 
the two modes |fc,+) and — ) at site k respectively. 
[T{k)] ~ ts{k)l + td{k)Tx is the 2x2 tunneling matrix 
that couples the two modes at site k to the two modes at 
site fc-|-l, Tx is the cc-Pauli matrix, and tg > {td > 0) de- 
note the tunneling amplitude for processes that preserve 
(flip) the local degree of freedom (Fig. [T]). We choose a 
parity-symmetric, real tunneling function t{k) = t{N—k) 
to ensure that Hq commutes with the combined PT op- 
erator pij. Note that we have chosen tunneling matrix 
so that Eq.([I]) is invariant under the exchange of pseu- 
dospin labels (T O —a. For a pair of balanced loss or gain 
impurities at mirror-symmetric locations, the potential is 
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FIG. 1. (color online) Tight binding lattice with a local degree 
of freedom represented by pseudospin cr = ±1. The tunneling 
amplitudes tsik) and tdik) denote processes that preserve or 
change the pseudospin respectively, while the tunneling from 
site k to fc + 1. Similarly, ±173 and ±17^ are on-site, non- 
Hermitian potentials that represent gain (white site) or loss 
(red site). This system is PT-symmetric irrespective of the 
time-reversal properties of the pseudospin. 

given by 

V = a'^-^iVam ~ ira^ (2) 

where fh = N + \ — m, i\r] — 175! + i^d^x denotes the 
non-Hermitian gain matrix at site m, and < 7^ < 7^ de- 
note the gain amplitudes for mode preserving and mode 
exchanging processes. The potential V is also invariant 
under the exchange of pseudospin labels, and is VT sym- 
metric irrespective of the time-reversal properties of the 
pseudospin. 

The eigenvalue difference equation obeyed by a two- 
component eigenfunction ^(k) — {fk,gk)'^ with energy e 
is given by 

- T{k - l)-^{k - 1) - T{k)m{k + 1) + 

{5k,^-5k,^)iT^{k) = e^{k). (3) 

where k = 1, • • • , iV. We note that open boundary con- 
ditions are implemented by assigning T(0) = = T{N) 
whereas periodic boundary conditions imply T(0) — 
T{N) 7^ 0. Using the symmetric and antisymmetric ba- 
sis that diagonalizes the tunneling matrix T{k) at every 
site, it is straightforward to obtain the following decou- 
pled equations, 

- (if/f+i + tLifti) + n''fi{Sk,m - 4,,^.) ,(4) 
-itkfk+i + ttjti) + *7^/.^(4,™ - 4,.^) =e/,^(5) 

Here t^^"^^ — [ts{k) ± td{k)] are the symmetric and anti- 
symmetric combinations of the tunneling rates, 7'^("^) — 
(7s ± 7d), and /f = {fk ± 9k) are the eigenfunction 
components in the symmetric-antisymmetric basis. 

Eqs.(|4|-([5| show that the T'T-symmetric Hamiltonian 
is a direct sum of Hamiltonians for two lattices with no 
internal structure: H — Hq + V — Hg Hj^ where Hs 
is the T'T-symmetric Hamiltonian with tunneling pro- 
file tf and a pair of non-Hermitian impurities at mirror- 
symmetric locations {m,rh) with strength 7'^, and Ha 
is obtained correspondingly. We emphasize that this de- 
composition into uncoupled problems is valid for arbi- 




FIG. 2. (color online) Tight binding ring with a local degree 
of freedom represented by pseudospin a = ±1 with a pair 
of "PT-symmetric impurities ±17 located at arbitrary sites 
(m,m). The tunneling matrix T is constant, but different, 
along the two paths that go from the gain site m (white site) 
to the loss site rh (red site). This system will map onto two, 
uncoupled rings with constant tunnelings = (tos ± Um) 

and if — {tbs ± tbd) along the two, outer and in-between, 
paths from site m to site m. 

trary, position-dependent tunneling profiles ts{k), mode- 
mixing amplitudes td{k), open or periodic boundary con- 
ditions, and arbitrary loss or gain strengths, as long as 
the underlying Hamiltonian is invariant under the ex- 
change pseudospin indices. 

Specific Cases and Numerical Results: When there is no 
mixing between the two pseudospin states, td ^ = 7d, 
the problem is trivial. In general, the T^T-symmetric 
threshold for Hq is equal to the smaller of the corre- 
sponding thresholds for Hs and Ha- 

When jd = 0, the loss (or gain) potential couples 
maximally to the pseudospin eigenmodes cr = ±1, and 
not to a linear combination of them. In this case, if 
the tunneling is constant, the T^T-symmetric phase di- 
agram ^pT{in) is given by a U-shaped curve, obtained in 
Ref. [TH], with the maximum value 7pt = [ts — td). For 
parity-symmetric, non-constant tunneling profiles, the 
appreciably strong VT- symmetric threshold, obtained 
in Ref. [52], is now selectively suppressed by increasing 
the mode-mixing tunneling amplitude td{k). For a lat- 
tice with periodic boundary conditions, we consider the 
model with tunneling matrices Tq = ipsl + todTx and 
Tb = tbsl + thdTx that are constant along each of the two 
paths that connect the gain site to the loss site (Fig.[2|. 
It then follows that the PT-symmetric threshold is inde- 
pendent of the distance between the loss and gain sites, 
as discussed in Ref. [2D], and is given by the smaller of 
the two combinations, {Iq — t^) and (<q — i^). Thus, 
a T^T-symmetric ring with a local degree of freedom of- 
fers significant threshold tunability independent of the 
distance between the loss and gain impurities. 
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FIG. 3. (color online) Typical PT-symmetric phase diagram for even (A'' = 40, left panel) and odd (A^ = 41 right panel) lattices 
with a local degree of freedom a. The vertical axis shows the PT-threshold "fprifJ.), measured in units of ts and the horizontal 
axis shows the fractional position ^ = m/N of the impurity. The tunneling is given by ts > td > 0, and impurity at site m acts 
as a gain for a = +1 and a loss for a — —1. When the mode-mixing td = (blue circles), decoupled pseudospins lead to earlier 
results [18]. As td increases (red squares and black stars), the PT-symmetric threshold 7pt(m) generally decreases. 



When 7d 7^ 0, the analysis carried out here predicts 
bounds on the gain matrix, given by (7s + 7d) < {tg + td) 
and (7s — 7^) < (is — td), however, these bounds do not 
determine the individual thresholds for 7s and 7^. In 
the extreme case of 7s = (meaning the gain poten- 
tial only couples to the symmetric combination), wc find 
that 7^ = 0, Ha is a purely Hermitian Hamiltonian and 
therefore, the PT-symmetric threshold is solely deter- 
mined by the Hamiltonian Hs- Note that, in general, a 
direct-sum decomposition of the Hamiltonian H is pos- 
sible if and only if the tunneling matrix T{k) at every 
site k and the non-Hermitian potential matrix iV can be 
simultaneously diagonalized. 

Finally we consider the case where the full Hamiltonian 
H cannot be decomposed into two non-interacting pieces. 
Generically, for an open lattice or a ring with constant 
tunneling matrix T and a single pair of gain/loss matrix 
iV , wave function matching approach [171 118) leads to a 
characteristic equation for eigenvalues of H that results 
from the determinant of a 6x 6 matrix. It is, thus, of little 
analytical value to calculate the T'T-symmetric threshold 
lPT{fn) and instead, we obtain the T'T-symmetric phase 
diagram numerically. We restrict ourselves to the sim- 
plest case of a constant-tunneling Hamiltonian Hq and 
an impurity potential matrix iT = ijsTz where Tz is the 
2;-Pauli matrix. In contrast to the previous cases, where 
the losses or gains for both modes occurred in the same 
waveguide, this non-Hcrmitian potential represents gain 
for one mode, a = +1, and loss for the other, cr = —1, at 
site m. 

Figure [3] shows the numerically obtained results for 



the threshold ^pT{l^)/ts as a function of the fractional 
location /i = m/N of the first impurity for three differ- 
ent values of mode-mixing tunneling td/ts = {0, 0.4, 0.7}. 
The left-hand panel shows the results for an even lat- 
tice with N = 40. When td = (solid blue circles), the 
two degrees of freedom are uncoupled and PT-symmetric 
phase diagram is identical to that for an open lattice 
with no internal degree of freedom [18]. As td/tg in- 
creases (solid red squares and black stars), generically, 
we find that the critical "fpTifJ-) is non-monotonically 
suppressed for different values of impurity locations /i. 
The right-hand panel shows corresponding results for 
an odd lattice with = 41. When td = (solid 
blue circles), the threshold impurity strength is given by 
ipr/ts = y/TTT/N « 1.012 when to = 1 :17], and 
therefore, does not appear in the figure. Once again, 
when td increases, the T^T-symmetric phase is (mostly) 
suppressed in a non-monotonic way. These results sug- 
gest that T^T-symmetry breaking in such systems shows 
a rich behavior that cannot be described with a simple 
analytical model. 

Discussion: In this paper, we have introduced VT- 
symmetric lattices with a local, two-state, quantum de- 
gree of freedom. By imposing invariance requirements 
on the Hermitian tunneling term, and PT-symmetric 
potential term that represents spatially separated gain 
and loss impurities, we have shown that a broad class 
of such lattice systems can be expressed as the direct 
sum of two, uncoupled, T'T-symmetric systems. In such 
cases, we have predicted that T^T-symmetric threshold 
can be tuned by mode-mixing tunneling amplitude. Since 
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the mapping is exact, all signatures of T-'T'-symmetry 
breaking, such as the ubiquitous, maximal chirality at 
■pT-symmetry breaking threshold ^U\, the even-odd ef- 
fect [35], tunable amplification f5S], etc. will be applica- 
ble in these cases as well. 

Since we have used the mode polarization as an ex- 
ample of the local degree of freedom, a microscopic cal- 
culation of the mode structure and the overlap between 



modes in adjacent waveguides is necessary to obtain typi- 
cal tunneling matrix elements. Similarly a detailed study 
of the selection rules for different polarizations will be 
necessary to characterize the relative strengths of ele- 
ments of the gain matrix 7^ and 7^. 
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